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A field-theoretical description of the behavior of compressible Ising systems with long-range inter- 
actions is presented. The description is performed in the two-loop approximation in three dimensions 
with the use of the PadBorel resummation technique. The renormalization group equations are ana- 
lyzed, and the fixed points that determine the critical behavior of the system are found. It is shown 
that the effect of elastic deformations on a system with a long-range interaction causes changes in 
its critical, as well as multicritical, behavior. 
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In compressible systems, the relation of the order pa- 
rameter to elastic deformations plays an important role. 
Earlier [1] it was shown that, in the case of an elasti- 
cally isotropic body, the critical behavior of a compress- 
ible system with a quadratic striction is unstable with 
respect to the relation of the order parameter to acous- 
tic modes, and a first-order phase transition close to a 
second-order one is realized. However, the conclusions 
formulated in the cited paper [1] hold only for low pres- 
sures. As was shown later [2], at high pressures begin- 
ning from a certain tricritical point Pt , the deformations 
induced by the external pressure affect the system to a 
greater extent and lead to a change of sign of the effective 
interaction constant for the order parameter fluctuations 
and, as a consequence, to a change in the order of the 
phase transition. In this case, according to [2], a homoge- 
neous compressible system is characterized by two types 
of tricritical behavior with a fourth-order critical point 
formed as the point of intersection of the two tricritical 
curves. Calculations performed in terms of the two-loop 
approximation [3] confirmed the presence of two types of 
tricritical behavior for Ising systems and provided values 
of the tricritical indices. 

In structural phase transitions that occur in the ab- 
sence of the piezoelectric effect, in the paraphase the elas- 
tic strains play the role of a secondary order parameter 
whose fluctuations are not critical in most cases. 

The effect of the long-range interaction described by 
the power law l/r~^~"' at long distances was studied 
analytically in terms of the e-expansion [46] and numer- 
ically by the Monte Carlo method [79] in two and three 
dimensions. It was found that the long-range interaction 
considerably affects the critical behavior of Ising systems 
for the parameter values a < 2. A recent study carried 
out for a three-dimensional space in the two-loop approx- 
imation [10] confirmed the prediction of the e-expansion 
for systems with long-range interactions. 

This paper describes the critical and tricritical behav- 
ior of three-dimensional compressible systems by taking 



into account the effect of the long-range interaction with 
different values of the parameter a. 

For a homogeneous Ising-like model with elastic defor- 
mations and a long-range interaction, the Hamiltonian 
can be represented in the form 



Ho = I d^xi^in + v-)sixr + +^isixrr] 
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a, (3=1 



(1) 



where S{x) is the scalar order parameter, uq is a positive 
constant, tq ~ |T — Td, Tj, is the phase transition temper- 
ature, Uaf3 is the strain tensor, ai and 02 are the elastic 
constants of the crystal, and aj, is the quadratic striction 
parameter. Let us change to the Fourier transforms of 
the variables in Eq. (1) and perform integration with re- 
spect to the components depending on the nonfluctuating 
variables, which do not interact with the order param- 
eter S{x) Then, introducing, for convenience, the new 
3 

variable y{x) = ^ Uaa{x), we obtain the Hamiltonian 
of the system in the form 



+ ^ J d^ qiSqlSq2Sq3S-ql-q2-q3 + 
+03 / d^qyqlSq2S-ql-q^, 



,(0) 

/ d'^qSqS-q 



1 f lai°^ 2 

+ 2^1 / d qyqV^q + -—Vq 



(2) 
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In Eq. (2), the components describing the homoge- 
neous strains are separated. According to [1], such a sep- 
aration is necessary, because the inhomogeneous strains 
Uq are responsible for the acoustic phonon exchange and 
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lead to long-range interactions, which arc absent in the 
case of homogeneous strains. For the system under study, 
let us determine the effective Hamiltonian that depends 
only on the strongly fluctuating order parameter S 



eM-H[S]} = B J exp{-HR[S,y]}Y[dyg 



(3) 



If the experiment is performed at constant volume, yo is 
a constant and the integration in Eq. (3) is performed 
with respect to only the inhomogcncous strains, while 
the homogeneous strains do not contribute to the effec- 
tive Hamiltonian. If the experiment occurs at constant 
pressure, the term Pfl is added to the Hamiltonian, the 
volume is represented in terms of the strain tensor com- 
ponents as 



UaaU0i3 + 0{U^)\ (4) 



and the integration in Eq. (3) is performed with respect 
to the homogeneous strains as well. According to [11], the 
inclusion of quadratic terms in Eq. (4) may be important 
when dealing with high pressures and crystals with large 
striction effects. The neglect of these quadratic terms re- 
stricts the applicability of the results obtained by Larkin 
and Pikin [1] to the case of low pressures. Thus, the 
Hamiltonian has the form 

+ ^(2:0 -W'o) j d^qid'^q2SqiS-qiSq2S-q2, (5) 

zo = al/{Aas), wo = a^^^^/{4a^^\ 

The effective interaction parameter vq = uq — Z()/2 that 
appears in the Hamiltonian due to striction effects, which 
are determined by the parameter zq, can take not only 
positive but also negative values. As a result, the Hamil- 
tonian describes both first-order and secondorder phase 
transitions. At vq = 0, the system exhibits a tricriti- 
cal behavior. In its turn, the effective interaction de- 
termined in Hamiltonian (5) by the parameter difference 
^0^0 may cause a second-order phase transition in the 
system when zqWq > and a first-order phase transition 
when zqWq < 0. This form of the effective Hamiltonian 
suggests that a higher order critical point can be realized 
in the system as the point of intersection of the tricritical 
curves when the conditions vq = and Zq = Wq are si- 
multaneously satisfied [2]. It should be noted that, with 
the tricritical condition zq = wq, Hamiltonian (5) of the 
model under consideration is isomorphic with the Hamil- 
tonian of a rigid homogeneous system. 

In the framework of the field-theoretical approach [12], 
the asymptotic critical behavior and the structure of the 
phase diagrams in the fluctuation region are determined 
by the CallanSymanzik renormalization group equation 



for the vertex parts of the irreducible Green functions. 
To calculate the /?— and 7-functions as the functions 
involved in the CallanSymanzik equation for renormal- 
ized interaction vertices u, ai, and a[o), or complex 



vertices z — af/4a3, 



q(o)2^^^(o)^ and v = ul2z, 
which are more convenient for the determination of crit- 
ical and tricritical behavior, a standard method based 
on the Feynman diagram technique and on the renor- 
malization procedure [13] was used with the propagator 
G(fc) — 1/t + jfej*^. As a result, the following expressions 
were obtained for the /3— and 7- functions in the two-loop 
approximation: 



/3„ = -{2a-D)v 1 - 36vJo + 1728(2Ji - - -G^v^ 



I3u 



^{2a-D)z 



1 - 24:vJo - 2zJo 



+ 576(2 Ji - - -G)v' 



-{2a-D)w 



1 - 24u Jo - AzJo + 2wJo 



+ 576(2 Ji - - -G)t;2 



It 



(2(7 - D) 



- 12i; Jo - 2z Jo -h 2w;Jo 



+ 288(2Ji- J2--G)t;2 

7y = (2£7-£>)192G^;^ 

f d^qd^p 

Ji = 

Jo = 
G = 



(1 + \q\-)^l + |p1-)(1 + +p2 + 2m1-/2) ' 
f d^q 

J a + w 



m 



o■^2 ' 



d 



d^grfXl + \q^ + k^ + 2fcg|'^) 

(1 + br )"'(! + I?' +p'+2pgr/')"' 



(7\-l 



Redefining the effective interaction vertices as 



Vi=V-Jo, V2 = Z-Jo, V3=W-Jo. (6) 



we arrive at the following expressions for the (3— and 7- 
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functions: 

/3i 



{2a-D)vi l-36t;i 

2 
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(7) 



+1728 



(2Ji 



p2 = -(2a-D)v2 l-24?;i-2w; 



+576(2Ji - 1 - -G)vl 



-{2a-D)v3 



1 - 2Avi - Av2 + 2v3 



+576(2Ji - 1 - -G)vl 



7t = (2fT - £)) - 12?;i - 2w2 + 2?;3 



+288(2Ji 



7y = (2c7 - Z))192Gt;i. 

The redefining proeedure makes sense for a a < D/2. 
In this case, Jq, J\, and G become divergent functions. 
Introducing the cutoff parameter A and considering the 
ratios Ji/Jq, G/Jq in the limit of A ^ oo, we obtain 
finite expressions. 

The values of the integrals were determined numeri- 
cally. For the case a a < i?/2, a sequence of the values 
of Ji/Jq and G/Jq was constructed for different A and 
then approximated to infinity. 

It is well known that perturbative series expansions 
are asymptotic and the interaction vertices of the order 
parameter fluctuations in the fluctuation region are suf- 
ficiently large to directly apply Eqs. (7). Therefore, to 
extract the necessary physical information from the ex- 
pressions derived above, the PadeBorcl method general- 
ized to the three-parameter case was used. The corre- 
sponding direct and inverse Borel transformations have 
the form 



J e *F{vit,V2t,vst,V4t)dt, 





F{vi,V2,V3) 



2 <^3 



For an analytical continuation of the Borel transform of 
a function, a series in an auxiliary variable is introduced: 



F{vuV2,V3,e)=Y,()'' Yl 



'-n,»2,»3 „tl„t2„,»3 



fe=0 



(8) 

and the [L/M] Fade approximation is applied to this se- 
ries at the point ^ = 1. This approach was proposed 
and tested in [14] in describing the critical behavior of 
systems characterized by several vertices of interaction 
of the order parameter fluctuations. The property [14] 
of the system retaining its symmetry under the Fade ap- 
proximants in the variable is essential in the description 
of multivertex models. 



In the two-loop approximation, the functions were cal- 
culated using the [2/1] approximant. The character of 
the critical behavior is determined by the existence of a 
stable fixed point satisfying the set of equations 



I3i{vt,v*2,v*s) = (i = 1,2,3). 



(9) 



The requirement that the fixed point be stable is reduced 
to the condition that the eigenvalues bi of the matrix 



dpi{vt,V*2,V*3) 

dvj 



(10) 



lie in the half-plane of the right-hand complex. The fixed 
point with v* = 0, which corresponds to the tricritical 
behavior, is a saddle point and must be stable in the 
directions determined by the variables z and w and un- 
stable in the direction determined by the variable v. The 
stabilization of the tricritical fixed point in the direction 
determined by the variable v is achieved by taking into 
account the sixth-order terms with respect to the order 
parameter fluctuations in the effective Hamiltonian of the 
model. The fixed point with z* = w* , which corresponds 
to the tricritical behavior of the second type, is also a sad- 
dle point and must be stable in the directions determined 
by the variables v and z and unstable in the direction de- 
termined by the variable w. Its stabilization is possible 
at the expense of the anharmonic effects. 

The resulting set of the resummed (3 functions contains 
a wide variety of fixed points. The table specifies the 
fixed points that are of most interest for describing the 
critical and tricritical behavior and lie in the physical 
region of vertex values with v, z, w > 0. The table 
also shows the eigenvalues of the stability matrix for the 
corresponding fixed points and the critical indices v and 
V- 

The analysis of the values and stability of the critical 
points suggests the following conclusions. Qualitatively, 
the critical phenomena seem to be identical for any value 
of the long-range interaction parameter a. The critical 
behavior of incompressible systems is unstable with re- 
spect to the deformation degrees of freedom (points 1). 
The stable point proves to be the one at a constant strain 
(points 2). Fixed points 3 describe the first type of tri- 
critical behavior of compressible systems, which occurs 
at constant pressure. Fixed points 4 are tricritical for 
systems studied at constant volume. Points 5 are fourth- 
order critical points at which two tricritical curves inter- 
sect. 

For the tricritical behavior of the first type (points 3), 
Hamiltonian (5) is isomorphic with the Hamiltonian of 
an incompressible homogeneous model and, hence, the 
critical indices also coincide with those of the incom- 
pressible model. The tricritical behavior of the second 
type (points 4) corresponds to the critical behavior of 
a spherical model and is determined by the correspond- 
ing indices. The fourth-order fixed points (points 5) are 
characterized by the field-average values of the critical 
indices. 
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The large values of the effective vertices z and w in 
comparison with the systems with short-range interac- 
tions [3] arc caused by the fact that the mechanism gov- 
erning the effect of elastic deformations on the critical 
phenomena is related to the dependence of the interac- 
tion integral in the Ising model on the distance between 
the lattice sites. 

The study described above revealed the considerable 



effect of elastic deformations on the critical behavior of 
systems with a long-range interaction. This effect mani- 
fests itself as a change in the values of the critical indices 
for Ising systems along with the appearance of multicrit- 
ical points in the phase diagrams of the substances. 

The work is supported by Russian Foundation for Ba- 
sic Research N 04-02-16002. 
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